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Abstract. 

We present constraints on the values of Cl m , n, erg, obtained from measurements of the X- 
ray luminosity function of galaxy clusters as compiled in EMSS, RDCS and BCS galaxy 
cluster samples. The values obtained fl„, = 0.25^'ij, n = -l^'ov an ^ 0-8 = 0.8^q , are 
in agreement with WMAP data. We then put constraints on the quoted parameters and 
the dark-energy equation-of-state parameter, w, using Chandra observations of 37 clus- 
ters with < z >=0.55 derived from 400 deg 2 ROSAT serendipitous survey and 49 bright- 
est z — 0.05 clusters detected in the All-Sky Survey. In the case of ACDM model, we 
get n m = 0.25^o 8 and cr 8 = 0.75^^, while in the case of the wCDM model, we get 

£2 m = 0.28!°|°, w = -1.25!°3°, and °~z = °- 8 -oo9- ° ur resu l ts are consistent with those 
from recent analyses of type la supernovae, cosmic microwave background anisotropies, 
the X-ray gas mass fraction of relaxed galaxy clusters, baryon acoustic oscillations and 
cosmic shear. The improvement in data quality from EMSS, RDCS, and BCS to Chandra 
observations leads to an improvement in the constraints even if not of the same entity of the 
improvement in data. 
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1. Introduction 

A cluster of galaxies is a large collection of 
galaxies held together by their mutual gravi- 
tational attraction. The galaxy clusters occu- 
pying the high-mass tail of this population 
provide a powerful and relatively clean tool 
for cosmology, since their growth is predomi- 
nantly determined by linear gravitational pro- 
cesses, their number density is exponentially 
sensitive to the size of the perturbations, and 
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hence can strongly constrain it. Starting with 
1990's, they have consistently indicated low 
values of Q„, (both from the baryonic frac- 
tion arguments (White et al. 1993) and mea- 
surements of the evolution in the cluster num- 
ber density (Eke et al. 1998; Borgani et al. 
200 1 )) and low values of crE (H enry & Arnaud 
1991; Reiprich & Boringer 2002; Schuecker 
et al. 2003) -a result since then confirmed by 

1 cr 8 is the amplitude of the mass density fluc- 
tuation power spectrum over spheres of radius 
8/i~'Mpc, and M 8 is the mean mass within these 
spheres 
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cosmic microwave background (CMB) studies, 
cosmic shear, and other experiments (Spergel 
et al. 2007; Fu et al. 2008). For precision's 
sake, cluster surveys in the local universe are 
particularly useful for constraining a combina- 
tion of the matter density parameter Q,„ and 
erg. Surveys that probe the cluster population 
at higher redshift are sensitive to the growth of 
density fluctuations, allowing one to break the 
£2 m -crg degeneracy that arises from local clus- 
ter abundance constraints (Eke, Cole & Frenk 
1996; Bahcall & Fang 1998). Recently, X-ray 
study of the evolution of the cluster mass func- 
tion at z — 0-0.8 have convincingly demon- 
strated that the growth of cosmic structure has 
slowed down at z < 1 due to the effects 
of dark energy, and these measurements have 
been used to improve the determination of the 
equation of state parameter (Vikhlinin et al. 
2008). 

The recipe of constraining cosmological 
parameters by means of clusters is composed 
of three ingredients: 1) The predicted mass 
function of clusters, n(M,z), as a function 
of cosmological parameters (erg, Q. M , wo, w a , 
etcfl The mass function can be calculated an- 
alytically (Press & Schechter 1974; Sheth & 
Tormen 2002; Del Popolo 2006) or fitting the 
results of numerical simulations (Jenkins et al. 
2001; Reed et al. 2003, Yahagi et al. 2004; 
Tinker et al. 2008). 

2) Sky surveys with well understood selection 
functions to find clusters, as well as a rela- 
tion linking cluster mass with an observable. 
A successful solution to the former require- 
ment has been to identify clusters by the X- 
ray emission produced by hot intracluster gas, 
notably using data from the ROSAT All-Sky 
Survey (RASS; Trumper 1993). The ROSAT 
Brightest Cluster Sample (BCS; Ebeling et al. 
1998, 2000) and ROSAT-ESO Flux Limited X- 
ray sample (REFLEX; Bohringer et al. 2004) 
together cover approximately two-thirds of the 
sky out to redshift z - 0.3 and contain more 
than 750 clusters. The Massive Cluster Survey 
(MACS; Ebeling et al. 2001, 2007)- which 



2 wo and w a represents the dark-energy equation- 
of-state parameter, w at z = and at generic redshift 
z, or expansion parameter a, respectively. 



contains 126 clusters and covers 55 per cent 
of the sky extends these data to z — 0.7. 
ROSAT 400 sq. degree survey (Vikhlinin et 
al 08), serendipitous cluster catalogue con- 
taining 266 groups/clusters. 3) A tight, well- 
determined scaling relation between survey 
observable (e.g. L x ) and mass, with minimal 
intrinsic scatter. 

The most straightforward mass-observable 
relation to complement these X-ray flux- 
limited surveys is the mass-X-ray luminosity 
relation. For sufficiently massive (hot) objects 
at the relevant redshifts, the conversion from 
X-ray flux to luminosity is approximately in- 
dependent of temperature, in which case the 
luminosities can be estimated directly from the 
survey flux and the selection function is iden- 
tical to the requirement of detection. A dis- 
advantage is that there is a large scatter in 
cluster luminosities at fixed mass; however, 
sufficient data allow this scatter to be quan- 
tified empirically. Alternative approaches use 
cluster temperature (Henry 2000; Seljak 2002; 
Pierpaoli et al. 2003; Henry 2004), gas fraction 
( Voevodkin & Vikhlinin 2004) or Yx parameter 
(Rravtsov et al. 2006) to achieve tighter mas- 
sobservable relations at the expense of reduc- 
ing the size of the samples available for analy- 
sis. The need to quantify the selection function 
in terms of both X-ray flux and a second ob- 
servable additionally complicates these efforts. 
In this paper, we use the observed X-ray lumi- 
nosity function to investigate two cosmologi- 
cal scenarios, assuming a spatially flat metric 
in both cases: the first includes dark energy in 
the form of a cosmological constant (ACDM); 
the second has dark energy with a constant 
equation-of-state parameter, w (wCDM). The 
theoretical background for this work is re- 
viewed in Section 2. Section 3 presents the re- 
sults and Section 4 the conclusions. 



2. Theory 

In this section, We derive an expression for the 
XLF (using now the mass function and M-T, L- 
T relations obtained in Del Popolo (2002), and 
Del Popolo el al. (2005) respectively) and then 
we get some constraints to Q m and n, by using 
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Fig. 1. (a) The 1 and 3 sigma credible regions 
of the marginalized posterior probability distri- 
bution P(Omega,n) (b) Marginalized posterior 
probability distribution P(£2,„). The dotted line 
is the 3 <r confidence interval (c) Marginalized 
posterior probability Distribution P(n). The 
dotted line is the 3 a confidence interval. 



EMSS, the ROSAT BCS, the RDCS (Borgani 
et al. 2001) samples. 

The mass function was calculated accord- 
ing to the model of Del Popolo (2006) in agree- 
ment with N-body simulations like those of 
Yahagi, Nagashima & Yoshii 2004). 
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Fig. 2. (a) The 1 , 2 and 3 sigma cred- 
ible regions of the marginalized pos- 
terior probability distribution P(n, erg). 
(b) Marginalized posterior probability 
distribution P(o"8). The dotted line is 
the 3 a confidence interval. 



In order to convert the mass function in 
a Luminosity function (XLF), first we con- 
vert mass into intra-cluster gas temperature, 
by means of the M - T x relation, and then 
the temperature is converted into X-ray lumi- 
nosity, by means of the L x - T x relation. The 
M - T x relation is obtained analytically us- 
ing the model of Del Popolo (2002), while 
the L x - T x relation is that obtained in Del 
Popolo, Hiotelis & Peharrubia (2005) based on 
an improvement of the Punctuated Equilibrium 
Model (PEM) of Cavaliere et al. (1997, 1998, 
1999) We performed the calculations assum- 
ing ACDM spectrum P(k) oc k"T 2 (k), where 
the transfer function T(k) is given by, BBKS 
(equation (G3)). The total number of X-ray 
clusters observed between luminosity and red- 
shift limits L\ < L < L u and z\ < z < z u , i.e. the 
cumulative luminosity function is given by: 

N(L h L u ,zi,z u ) = f ' f"A(L,z) 
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Fig. 3. (a) The 68% and 95% 
credible regions of the marginal- 
ized posterior probability distribu- 
tion P(o- 8 ,fl m ) for a ACDM model, 
(b) Same as panel (a) but for 
a wCDM model, (c) The 68% 
and 95% credible regions of the 
marginalized posterior probability 
distribution P(w,fi,„) for a wCDM 
model. 



dn{L\,£) 
d~L 



dLdV(z) 



(1) 



where A{L\ , z) is the area of the sky that an X- 
ray survey samples at redshift z as a function 
of luminosity L. The previous model, consists 



of the parameters: h, Q.bh 2 , Q c /z 2 0, cr 8 , n, w, 
A\. We use priors as in Mantz et al. (2007) for 
h, Clbh 2 , n, Ai, and we obtain constraints for 
Q m , n, <x 8 , and w comparing model and data 
through a Bayesisn inference analysis. In or- 
der to perform the quoted analysis, we need 
the posterior probability distribution for Q m , 
n, erg, namely P(£2 m ,n,<r 8 ) which is obtained 
by normalizing the product of the prior proba- 
bility distribution and the likelihood function 
(Gregory & Loredo 1992). In order to con- 
straint, Q m , n and cr 8 , we use the likelihood 
function L{Q. m , n, cr 8 ) (see Cash 1979, Reichart 
et al. 1999): 



L(£l m , n, cr 8 ) = Y\ P(Li, Zi\£l m , n, erg) 



(2) 



where P(L,, z,if2 m , n, cr 8 ) is the probability that 
the 2th X-ray cluster fits our model, given val- 
ues of Q m , n, cr 8 . For our model, this probabil- 
ity is given by (e.g., Cash 1979): 



P(L,z\Cl m ,n,crs) = 



A(L u z) 



dn c (L,z) dV(z) 



dL 



dz 



N{L u L u ,zi,z, u ) 



(3) 



In Eq. (O, N tot is the total number of X- 
ray clusters in the same region of the L — z 
plane as that over which N{L\,L u ,z\,z, a ) is 
defined. The posterior probability distribution 
for one of the parameters, e.g., P(Q. m ), is ob- 
tained by marginalizing the posterior proba- 
bility distribution for all the three parame- 
ters, P(£l m , n, cr 8 ), over the other parameters 
(Gregory & Loredo 1992). The credible re- 
gions (1, 2, 3 cr) are determined by integrat- 
ing the posterior probability distribution over 
the most probable region of its parameter space 
until 68,3%, 95.4% and 99.73% (respectively) 
of this distribution has been integrated. 

3. Results 

The constraints obtained from the compari- 
son of the model of the previous section with 
EMSS, RDCS, and BCS are shown in Fig. 1- 
Fig. 2. In Fig. la, we plot the 1 and 3 cr con- 
fidence regions of the marginalized posterior 



3 Q b and £l c are the baryon and CDM densities 
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probability distribution P(O m ,n) for a ACDM 
power spectrum, while in Fig. lb, lc, the 
marginalized posterior probability distribution 
P(Q m ), and the marginalized posterior proba- 
bility distribution P(n). The dotted line repre- 
sents the 3 cr confidence interval. Fig. 1 shows 
that Q m = 0.25+°;^ and n = -l^os- Similarly, 
in Fig. 2a, we plot the 1, 2 and 3 cr credible 
regions of the marginalized posterior probabil- 
ity distribution P(n, <x 8 ), while in Fig. 2b the 
marginalized posterior probability distribution 
P(cr 8 ). As in Fig. 1, the dotted line represents 
the 3 cr confidence interval. The constraint that 
we obtain for <x 8 is <x 8 = 0.8+5}. Using the 
same model, we have performed the same anal- 
ysis for the ACDM and wCDM models us- 
ing more recent data, namely Chandra obser- 
vations of 37 clusters with < z >=0.55 de- 
rived from 400 deg 2 ROS AT serendipitous sur- 
vey and 49 brightest z - 0.05 clusters detected 
in the All-Sky Survey. Fig. 3a plots the 68.3 
and 95.4 % confidence constraints on <x 8 and 
Q m for the ACDM model. The constraints from 
the cluster sample used are Q,„ = 0.25 ^} ( ' 8 
and cr 8 — 0.75^qqj. These constraints are in 
good agreement with recent independent re- 
sults from the CMB (Spergel et al. 2007) and 
cosmic shear, as measured in the 100 Square 
Degree survey (Benjamin et al. 2007) (Figure 
8) and CFHTLS Wide field (Fu et al. 2008). 
Our results are also in good overall agreement 
with previous findings based on the observed 
X-ray luminosity and temperature functions of 
clusters (e.g. Eke et al. 1998; Schuecker et 
al. 2003; Henry 2004), Our result on Q m is 
in excellent agreement with current constraints 
based on cluster / gas data (Allen et al. 2008 
and references therein) and the power spectrum 
of galaxies in the 2dF galaxy redshift survey 
(Cole et al. 2005) and Sloan Digital Sky Survey 
(SDSS) (Eisenstein et al. 2005; Tegmark et al. 
2006; Percival et al. 2007), as well as the com- 
bination of CMB data with a variety of external 
constraints (Spergel et al. 2007). 

In the case of the wCDM model, Fig. 3a 
shows the joint constraints on Q„, and <x 8 from 
the luminosity function data, while Fig. 3b dis- 
plays constraints on Q m and w. The marginal- 
ized results from the X-ray luminosity func- 
tion data are Q m = 0.28^ }° w = -1.25!"" 



and <x 8 — 0.8!|j|j(j. XLF results are consis- 
tent with each WMAP, SNIa data, and cluster 
/gas data, and with the cosmological-constant 
model (w = -1). The results for <x 8 presented 
here are somewhat different than those from 
previous work based on the BCS and REFLEX 
data. The magnitude of this discrepancy under- 
scores the need for an improved understand- 
ing of the observational biases resulting from 
asphericity, projection effects and hydrostatic 
disequilibrium. More advanced and compre- 
hensive simulations, calibrated by gravitational 
lensing studies, show considerable promise in 
this area. With regard to dark energy, the 
systematics due to the theoretical mass func- 
tion and the redshift evolution of the mass- 
observable relation are also important. The for- 
mer can be addressed with a large suite of cos- 
mological simulations, The latter necessitates 
a more rigorous study of galaxy cluster virial 
relations and their evolution 

4. Conclusions 

In this paper, we presented constraints on the 
values of Q. m , n, cr 8 , and w obtained from 
measurements of the X-ray luminosity func- 
tion of galaxy clusters as compiled in EMSS, 
RDCS and BCS galaxy cluster samples. The 
values obtained using EMSS, RDCS, and BCS 
data are £l m = 0.25^, and n = -1 
cr 8 = 0.8+°}, in agreement with WMAP 
data. Using Chandra observations of clusters 
from 400 deg 2 ROSAT serendipitous survey, 
we get Q m = 0.25!^ g and <r s = 0.75™*, 
for a ACDM model, while in the case of the 
wCDM model, we get £l m = 0.28 + °-}°, w = 
-1.25+°;^, and <x 8 = 0.8+°;°*, in agreement 
with analyses of type la supernovae, cosmic 
microwave background anisotropics, the X-ray 
gas mass fraction of relaxed galaxy clusters, 
baryon acoustic oscillations and cosmic shear. 
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G. S. BISNOVATYI-KOGAN: What about 
the acoustic oscillations in the mass function: 
can you reveal them? 

A. DEL POPOLO: Baryon Scoustic 
Oscillations (BAOs) are mostly viewed in 
the two-point redshift-space correlation func- 
tion and its Fourier transform (the spectrum). 
Eisenstein (2005) detected the acoustic peak 
in the SDSS LRG survey, while Cole et al. 
(2005) discovered them in the 2DF spectrum. 
BAOs detected in 2dFGRS and SDSS have 
been used to constrain cosmological models 
by Percival et al. (2007). 
J. BECKMAN: Am I correct in seeing that 
the value of n in your results is not consistent 
with the value 1 ? 

A. DEL POPOLO: The results of the model 
using older data (EMSS, RDCS, BCS) are 
consistent with the value 1 (n = -l^os)- 



